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ABSTRACT 

We observed the Seyfert 1 galaxy NGC 4151 on eleven occasions at 1-2 day intervals using the Berkeley spec- 
trometer during the ORFEUS-SPAS II mission in 1996 November. The mean spectrum covers 912-1220 A at ^0.3 
A resolution with a total exposure of 15,658 seconds. The mean flux at 1000 A was 4.7 x 10~ 13 erg cm" 2 s" 1 A" 1 . 
We identify the neutral hydrogen absorption with a number of components that correspond to the velocity dis- 
tribution of H I seen in our own Galaxy as well as features identified in the C IV A 1549 absorption profile 
by Weymann et al. The main component of neutral hydrogen in NGC 4151 has a total column density of 
log Nhi = 18.7 ± 1.5 cm" 2 for a Doppler parameter b = 250 ± 50 km s" 1 , and it covers 84 ± 6% of the source. 
This is consistent with previous results obtained with the Hopkins Ultraviolet Telescope. Other intrinsic far-UV 
absorption features are not resolved, but the C III* Al 176 absorption line has a significantly higher blueshift rela- 
tive to NGC 4151 than the C III A977 resonance line. This implies that the highest velocity region of the outflowing 
gas has the highest density. Variations in the equivalent width of the C III* A 1176 absorption line anticorrelate 
with continuum variations on timescales of days. For an ionization timescale of <1 day, we set an upper limit of 
25 pc on the distance of the absorbing gas from the central source. The O VI A 1034 and He II A 1085 emission 
lines also vary on timescales of 1-2 days, but their response to the continuum variations is complex. For some 
continuum variations they show no response, while for others the response is instantaneous to the limit of our 
sampling interval. 

Subject headings: galaxies: active — galaxies: individual (NGC 4151) — galaxies:nuclei — galaxies:Seyfert — 
ultraviolet: galaxies 



1. INTRODUCTION 

The Seyfert 1 galaxy NGC 4151 has been a fav< 
servers wishing to use variability as a probe of the 
ings of active galactic nuclei (AGN). NGC 415Tj 



i-h , variability was first noted in a rocket flight ( Hartig 



i )rite of ob- 
i nner work- 
ultraviolet 
19791). Nu- 



merous campaigns using the International Ultraviolet Explorer 
(IUE) have since monitored the continuum, emissic n lines and 
absorption lines (|Ulrich et al 



1984 



Bromage et 



Clavel et al. 1992| ; penshaw et al. 1996fc pdelson e t al 1996| ) 
IUE studies of the absorption lines ( promage et al . 985) noted 



a tendency for the equivalent widths of the high ioni nation lines 
(e.g., C IV) to correlate directly with variations in the continuum 
while low ionization lines such as C II were antico related. In 
contrast, monitoring of the C IV absorption line a c ecade later 
using the Goddard High Resolution Spectrograph i GHRS) on 
the Hubble Space Telescope (HST) at moderate s jectral res- 
olution (R^ 15,000) showed no variation in streng :h or shape 



over a baseline of several years (Weymann et al. 19 )J) 



Significant UV continuum and emission-line variability 
timescales as short as 1-2 days were seen in the 
monitoring of N GC 4151 with IUE in 1994 December 
shaw et al 



1996). Monitoring with the Hopkins 
Telescope (HUT) during the Astro-2 mission ove 
1820 A band at ~ 2-day intervals showed no change 
equivalent width, but revealed substantial variation 
column and in low-ionization ions such as C III anil 
were correlated with variations in the continuum 
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1996; pCriss et al. 1998[ ). With the low resolution (R-300) 
HUT data, however, we could not identify which of the several 
comj onents in the complex C IV absorption profile was associ- 
ated vith these variations. 

The HUT Astro- 1 and Astro-2 observations of NGC 4151 
led numerous absorption lines ov er a wide range o f ion- 
n stat es in the 912-1200 A band dKriss et al. 1992]; Kriss 
1995). Our observations with the Berkeley spectrometer 
ORFEUS-SPAS II were intended to make use of its higher 
resol ition (R~3300) relative to HUT to study the structure and 
time response of the absorption line systems in greater detail. 
We j resent here the mean spectrum and key features of the 
variability we saw during the ORFEUS-SPAS II mission. Full 
detai s of the individual observations will be described in a sub- 
sequ< :nt publication. 



2. OBSERVATIONS AND DATA REDUCTION 

obtained spectra of NGC 4151 at eleven epochs with the 
:ley spectrograph during the ORFEUS-SPAS II mission in 
November/December. Exposure times ranged from 600 
s, with a total exposure of 15,658 s. The general de- 
af the Berkeley spectrograph is discussed by Hurwitz & 
'er (1986, 1996), while calibration and performance for 
RFEUS-SPAS II mission are described by Hurwitz et al. 



The Berkeley spectrometer has a 2-D detector, so the 



'ORFEUS-SPAS II Guest Investigator 
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NGC 4151 data were obtained through a 26" diameter aper- 
ture simultaneously with airglow spectra through a larger aper- 
ture offset by 2.4'. The signal-to-noise ratio varies greatly with 
wavelength, and it is highest near the O VI and Lya emission 
lines. Due to a reflection in the Berkeley spectrometer, undis- 
persed Lya airglow produces a rapid increase in the background 
level below ~950 A. Correction for this and other instrumen- 
tal features are discussed in Hurwitz et al. (1997), and our 
data were extracted and calibrated using the prescription rec- 
ommended by the Berkeley spectrometer team. 

To maximize the signal-to-noise ratio (S/N), we binned these 
basic spectra into larger wavelength bins of 0.15 A, roughly 
half a resolution element ( fr-lurwitz et al. 1997a ). To correct 
the wavelength scale for slight inaccuracies in telescope point- 
ing, we shifted the object spectra to place the Galactic Si II 
Al 190 and Al 193 absorption features at a heliocentric velocity 
of 6 km s" 1 (the Galactic H I velocity determined by Murphy 
1997). From a comparison among all the observed Galactic 
absorption features in the spectrum, we estimate that our wave- 
length scale is accurate to 40 km s" 1 (lc). 

The contemporaneous airglow spectra were extracted in ex- 
actly the same manner as the target data. Airglow emission 
lines were fitted using Gaussian components. Comparison of 
the line widths of isolated emission lines through both the air- 
glow aperture and the smaller target aperture provided the cor- 
rection factor needed to generate a model airglow spectrum tai- 
lored to each target spectrum. Our template airglow spectra 
were then corrected for slight wavelength offsets and scaled 
to match the flux observed through the target aperture by lin- 
ear shifting and scaling to the Ly/3 line in the object spectrum. 
The mean airglow-subtracted, exposure-weighted spectrum is 
shown in Figure 1 . 

The mean spectrum shows that NGC 4151 had an intensity 
intermediate between the observations made during the Astro- 1 
and Astro-2 missions. We use the mean spectrum to determine 
the basic parameters of both continuum and line components 
(see Tables 1 and 2). These and subsequent fits were performed 
using \ 2 min imization in the spectral fitting program spec fit 



Table 1 

Absorption lines in the mean spectrum of NGC 4151 



(Kriss 1994). We model the continuum as a power law in f\. 
We correct for foreground extinction using a Clayton, Cardelli, 
& M athis (1989) curve assuming E(B-V) = 0.04 ( Kriss et al. 
1995) and R v = 3.l. 

The best-fit, extinction-corrected continuum is /a = 8.20 x 
10" 13 (A/T000)" 137 . For the Lya, He II A1085 and O VI A1036 
emission lines we assume the shape found by Kriss et al. 
(1992). Absorption features were modeled using components 
that are either Gaussians in equivalent width or optical depth. 
The features identified in the mean spectrum and their parame- 
ters are listed in Table 1 . We find no evidence for a change in 
power-law index among the individual observations. To fit the 
variations in the individual emission and absorption lines, we 
fix the wavelength and line width of each component and fit the 
intensity or equivalent width in every individual exposure. 

3. ABSORPTION FEATURES AND THEIR VARIABILITY 

3.1. Neutral Hydrogen Absorption 

Measuring the intrinsic neutral hydrogen column in 
NGC 4151 requires good S/N in the high-order Lyman lines 
since the low-order lines are heavily saturated. A S/N of ~ 5 
per bin at the Lyman limit required further binning of the the 
mean spectrum to 0.6 A pixel" 1 . This is still far better than the 
~3 A resolution of HUT. Unfortunately, the high background 
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a The Galactic Si II A 
align the wavelength 
b N: lines intrinsic to 



1 190 and Al 193 lines (with vq = 60 km s 1 ) were used to 
scale. Velocities are accurate to ~ 40 km s" 1 (1<t). 
NGC 4151; G: Galactic absorption lines. 



at short wavelengths prevents a good measure of the hydrogen 
column in the individual observations. No significant varia- 
tions in the strong Lya and Ly/3 absorption lines are seen, pre- 
sumably because of their high optical depth. The higher spec- 
tral resolution of the Berkeley spectrometer compared to HUT, 
however, enables us to deblend the hydrogen absorption in the 
mean spectrum and separate the complex assortment of Galac- 
tic and intrinsic components that are present. 

The outflowing H I absorption in NGC 4151 brackets strong 
absorption by neutral hydrogen in our own Galaxy. Galactic 
H I 21 -cm emission observed toward NGC 4151 shows multi- 
ple components with a total column density of ~ 2 x 10 20 cm" 2 
at a mean heliocentric velocity of 60 km s" 1 (M urphy et al. 
1996; Murphy 1997| ). This complex structure is reflected in our 
ORFEUS-SPAS II data by the large line widths observed for 
Galactic absorption features. 

To model the Galactic H I absorption, we first fit the absorb- 
ing column, location and line width of the four main compo- 
nents in Murphy's 21 -cm profile. We then generated an H I 
absorption line template smoothed to the resolution of our spec- 
trum. We find that although the total absorbing column is 
2 x 10 2() cm" 2 and the separation of the components is less than 
our instrumental resolution, the resulting absorption line profile 
is optically thick but not black, consistent with our observations 
of the high-order Lyman lines. 

We similarly determined the location and line width of 
two narrow low-ionization components seen in GHRS spec- 
tra (components A and F) which are identified with outflow- 
ing and halo material in NGC 4151, respectively ( Weymann et 
al. 1997). We fit these components in the mean spectrum and 
included a broad component to represent higher ionization ma- 
terial (an unresolved blend of components C, D and E). The 
several components and their best fit parameters are summa- 
rized in Table 2. The reduced \ 2 of the fit is quite large (187.6), 
but this is largely due to an underestimate of the true errors in 
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FIG. 1. — Observed mean airglow-subtracted spectrum of NGC 4151 obtained with the Berkeley spectrometer on the ORFEUS-SPAS II mission in Novem- 
ber/December 1996. 



the binned data. Near O VI in the re-binned spectrum the formal 
S/N based on photon statistics is ~85, but the fiat field is only 
accurate to ~ 5%. 

Our best fit to the broad absorption gives log(NHi) =19.7 
and b = 200 km s" 1 covering ~ 84% of the source. This broad 
component is co nsistent with that seen in the HUT A stro- 1 and 
Astro-2 spectra ( |Kriss et al. 1992| ; [Kriss et al. 1995| ). It has a 
width and partial covering similar to that seen during Astro-2, 
and a total column intermediate between the Astro- 1 and Astro- 
2 values, consistent with the intermediate value of the contin- 
uum flux seen in these new observations. Component 'F', at the 
systemic velocity of NGC 4151 and probably associated with 
its ISM or halo, has a column density an order of magnitude 
lower than the optically thick component suggested by Kriss et 
al. (1995), but it is still optically thick. Given the lack of broad, 
black Mg II absorption at this velocity in the GHRS spectrum 
(Weymann et al. 1997), it is likely that the actual column den- 



sity is at the lower end of the range permitted by our errors. 

3.2. Other Absorption Features 

Correcting for the 90 km s" 1 resolution of the Berkeley 
spectrometer, the unresolved Galactic absorption features have 
FWHM ~ 100 kms -1 (see Table 1). This is roughly what is ex- 
pect ed from the 21 -cm H I profile towards NGC 4151 ( Murphy 
1997), and is consistent with Galactic lines in the GHRS spec- 
tra. Our identified Galactic absorption lines are similar to those 
reported for the 3C 273 sight line by Hurwitz et al. (1997b). 



Accurate measurements of O VI absorption intrinsic to 
NGC 4151 proved difficult due to the overlapping absorption 
and emission features. What is apparent, however, is that the 
dominant absorption in the O VI doubl et is similar to the ab- 



sorption features seen in the C IV line (Weymann et al. 1997 
Hutchings et al. 1997] ). It is optically thick, but the broad width 
(~ 1100 km s" 1 ) and the non-black line centers imply partial 
covering. We see no significant variations in the O VI ab- 
sorption equivalent width. More detailed analysis of the O VI 
troughs will be reported in a future publication. 

As described by Bromage et al. (1985) and Kriss et al. 
(1992), C III* A1176 can be used together with C III A977 to 
measure the density of the absorbing gas. In the Astro- 1 and 
Astro-2 HUT data, both lines were assumed to arise in the same 
absorbing region. The higher resolution ORFEUS-SPAS II data 
shows that the situation is not so simple — the mean velocity of 
the C III* Al 176 absorption is significantly higher than that of 
C III A977. C III A977 is less blueshifted, and it has a veloc- 
ity comparable to that of component D in the C IV profile of 
Weymann et al. (1997). C III* Al 176 is closer in velocity to 
component C, but the 850 km s" 1 width of the feature gives 
considerable overlap to all broad C IV components C, D, and 
E. If the gas outflowing from NGC 4151 is accelerating, this 
places the highest density gas in the outflow farthest from the 
source. 

In contrast to the results for Ly/3 and O VI absorption, we 
detect strong variations in the C III* Al 176 and S IV A 1073 ab- 



Table 2 
Hi Absorption in NGC 4151 



Component 


vq (kms ') 


log N(HI) (cnT 2 ) 


b(kms~') 


Covering fraction 


NGC 4151 'A' 


-540 


15.5 ±0.1 


45±0 


0.95 ±0.24 


Galactic 


60 


20.3 ±0 


75±0 


1.00±0 


NGC 4151 'C+D+E' 


451 


18.70± 1.5 


250 ±50 


0.84 ±0.06 


NGC 4151 'F' 


982 


19.61 ±0.5 


45 ±0 


0.73 ±0.80 
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sorption lines. (S IV A1062 is blended with Galactic absorp- 
tion.) The constancy of the Galactic features throughout the 
mission confirms the reality of these variations. As shown in 
Figure Q both C III* and S IV anti-correlate with the observed 
continuum variations. To within the 1-2 day sampling of our 
measurements, the continuum and line variations are instanta- 
neous with no significant lag. 

Following the formalism of Krolik & Kriss (1997), we can 
use the ionization and recombination timescales to set limits on 
the location of the absorbing material. For an observed ioniza- 
tion timescale r lo „, an ionizing flux density at the earth /} „, a 
mean photoionization cross section (a km), an ionization energy 
hvj, and a distance to the source D (15 Mpc for NGC 4151), 
they find that the material must lie at a radius 



r < 



fion \&h 



llVj 



1/2 



D. 



(1) 



The decrease in the C III* Al 176 equivalent width during the 
increase in the continuum flux observed from day 334 to 335 
implies that ionization is dominating in this interval. We thus 
have an upper limit of 1 day on f,- on . Usi ng a photoionizat ion 
cross section of (cr,„„) = 1.6 x 10" 18 cm 2 ( psterbrock 1989 ) at 
the ionization threshold for C III* (41.4 eV, or 300 A), and ex- 
trapolating our best-fit, extinction corrected continuum to 300 
A, we find r < 25 pc. This is comparable to the limit found us- 
ing the variation in the H I col umn seen during the H UT Astro- 
2 ob servations of NGC 4151 flKrolik & Kriss 1997|; Kriss et al. 
1998). 

4. EMISSION FEATURES 

The dominant emission features in our data are O VI A 1034, 
He II A1085, and the blue wing of Lya. The partial Lya line 
and the strength of the Lya airglow preclude any meaningful 
determination of the intrinsic line strength. Variations in the 
O VI and He II fluxes compared to the continuum are shown in 
Figure |j. The fractional variation in the He II flux (F var = 0.23) 
is comparable to that in the continuum flux (f ra( -(1000) = 0.20). 
In contrast, the variation seen in the O VI emission is much 
less: F var = 0.12. Our data set is not large enough to com- 
pute meaningful cross correlation functions, but qualitatively 
one can see that there is significant response in the emission 
lines on timescales of 2 days or less. The detailed response, 
however, is not the typical smeared and delayed variation seen 
in prior monitoring campaigns. While the continuum is rising 



in the first three observations, O VI and He II are constant or 
falling in flux. In contrast, both lines track well the dip and 
subsequent rise in continuum flux at days 334 to 336. 

We thank the Berkeley spectrometer science team and all 
those whose participation helped make the ORFEUS-SPAS II 
mission such a success. This work was supported by NASA 
LTSA grant number NAG 5-3255 to the Johns Hopkins Univer- 
sity. 
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FIG. 2. — Absorption line, continuum, and emission line variations with time 
for the eleven observations of NGC 4151 during the ORFEUS-SPAS II mis- 
sion. Note that the Galactic Si II A 1190 equivalent width in the top panel 
shows no evidence for variation. The continuum flux at 1000 A is in units 
of 10~ 13 erg cm" 2 s" 1 A" 1 . The He II A1085 and O VI A1034 emission line 
fluxes are in units of 10~ 13 erg cirT 2 s _1 . 
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